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Abstract 

The calculation of the 1 = 2 pion scattering length in quenched lattice 
QCD is revisited. The calculation is carried out with the Wilson fermion 
action employing Liischer's finite size scaling method at (3 = 5.9, 6.1, and 6.3 
corresponding to the range of lattice spacing a ~ 0.12 — 0.07 fm. We obtain 
in the continuum limit a^/m^ = —2.09(35) 1/GeV 2 , which is consistent with 
the prediction of chiral perturbation theory ao/m n = —2.265(51) 1/GeV 2 . 

PACS number(s): 12.38.Gc, 11.15.Ha 
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Lattice calculations of S-wave scattering lengths of the two-pion system are important 
step to understand dynamical effects of strong interactions. There are already a number of 
calculations for the 1 = 2 process with either Staggered [1,2] or Wilson fermion action [1,3]. 
While these calculations gave results that are in gross agreement with the prediction of 
chiral perturbation theory (CHPT) [4], they were made on coarse and small lattices. More 
importantly, the continuum extrapolation was not made. Aiming to improve on these points, 
we carried out a calculation of the 1 = 2 S-wave scattering length in quenched lattice QCD. 
A preliminary result was reported in Ref. [5], in which some disagreement with the CHPT 
prediction was mentioned. In the mean time Liu et. al. carried out a similar calculation with 
the improved gauge and the Wilson fermion actions on anisotropic lattices [6]. 

We employ the standard plaquette action for gluons and the Wilson action for quarks, 
and explore the parameter range m^/rrip ~ 0.55 — 0.80 for the chiral extrapolation and 
a ~ 0.07 — 0.12 fm for the continuum extrapolation. This is compared with the parameters 
of Liu et.al. which range m^/rrip ~ 0.7 — 0.9 and a s ~ 0.2 — 0.4 fm. Our calculations are 
made for parameters significantly closer to the chiral limit. In this short report we give the 
final result of our analysis. 

The numbers of configurations (lattice sizes) are 187 (16 3 x 64), 120 (24 3 x 64), and 
100 (32 3 x 80) for (3 = 5.9, 6.1, and 6.3, respectively. Quark propagators are solved with the 
Dirichlet boundary condition in the time direction and the periodic boundary condition in 
the space directions. The pion mass covers the range of 450 — 900MeV. The lattice constant 
is estimated from the p meson mass, which was obtained in our previous study [7], to be 
a' 1 = 1.64(2), 2.29(4), 3.02(5) (GeV) at (3 = 5.9, 6.1, and 6.3. Our calculations were carried 
out on the Fujitsu VPP500/80 supercomputer at KEK. 

The energy eigenvalue of a two-pion system in a finite periodic box L 3 is shifted by the 
finite-size effect. Liischer presented a relation between the energy shift AE and the S'-wave 
scattering length a , given by [8] 

—AE • = T + d • T 2 + C 2 • T 3 + 0(T 4 ) , (1) 

where T = a /(nL). The constants are C\ = —8.9136 and C 2 = 62.9205 computed from 
geometry of the lattice. Since T has a small value, typically ~ — 10~ 2 in our simulation, we 
can safely neglect the higher order terms 0(T 4 ). 

The energy shift AE can be obtained from the ratio R(t) = G(t)/D(t), where 

G(t) = (7r + (t)7r + (t)W-(ti)W-(t 2 )) 

D(t) = (n+(t)W-( tl )) (n + (t)W-(t 2 )) . (2) 

In order to enhance signals against the noise we use wall sources for 7r~, which are denoted 
by W~ in (2), by fixing gauge configurations to the Coulomb gauge. The two wall sources 
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are placed at different time slices t\ and t 2 to avoid contaminations from Fierz-rearranged 
terms in the two-pion state which would occur for ti = t 2 . We set t 2 — ti + 1 and ti = 8, 
10, 13 for = 5.9, 6.1, 6.3. 

An example of R(t) is plotted in Fig. 1 for f3 — 6.3 and k = 0.1513 corresponding to 
rrij, = 433 (4) Me V. We see a clear, almost linear fall-off as a function of t till t — 80 even 
for a small energy AE 20 MeV, showing that our wall sources work well for the two-pion 
state. 

The energy shift AE is obtained from the linear term in the expansion of R(t) : 

R(t) = Z ■ (l - AE ■ t + 0(r 2 )) , (3) 

where r = t — t 2 . The quadratic and higher order terms have no simple relations to AE due 
to effects from intermediate off-shell two-pion states [2] and quenching effects [9]. We first 
attempt to fit the data with the form 

(Sqr) Z ■ (l - AE ■ t + E' ■ r 2 ) . (4) 

We find that this fit (Sqr) is quite ill-determined, since the two terms correlate so strongly, 
resulting in unacceptably large errors in AE and E'. We then attempt to fit with 

(Exp) Z ■ exp(-A£ • r) , 

(Lin) Z ■ ( 1 - AE ■ r) . (5) 



These fitting forms give well-determined AE, while it may be contaminated by contributions 
from the second order term. We also include a fit of the form 

(Old) Z-AE-t (6) 

into our attempts for completeness, since this was used in our preliminary report [5]. Note, 
however, that this form is theoretically correct only when Z is close to unity. The results 
for AE (and E 1 in case (Sqr)) are given in Table I for (3 = 5.9, Table II for (3 = 6.1, and 
Table III for (3 = 6.3. We take the same fitting range for the four fits, t = 21 — 42 for (3 = 5.9, 
t = 25 - 50 for (3 = 6.1, and * = 27 - 62 for f3 = 6.3. The value of X 2 for each fitting is 
always small, and does not discriminate among fits. We do not consider case (Sqr) further 
because of very large errors, although the resulting central values for the energy shift are 
consistent with those from (Exp) and (Lin). The problem we must consider is whether we 
can remove contaminations of the second order term for AE from (Exp) and (Lin). 

Figure 2 shows ao/m w as a function of the pion mass obtained at each j3, with their nu- 
merical values tabulated in Table I, II, and III. We observe a large difference between (Exp) 
and (Lin), indicating that contributions from the 0(r 2 ) term are indeed non- negligible and 
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largely affect the determination of AE. The common in all figures of ao/m n versus m n is 
that the data show a behavior linear in ml. We then fit 



a /m n = A + B ■ m 2 



(7) 



to extract the value A in the chiral limit. From the view point of CHPT we may in principle 
have a term m 2 log(m 2 /A 2 ) added to (7). If we include this term with a free coefficient into 
the fit, however, the coefficients correlate too strongly that the fit is invalidated, producing a 
large error also for A. It is difficult to distinguish m 2 and m 2 log(m 2 ) within the range of m 2 
that concerns us and the limited statistics. Since we do not see any significant curvature in 
the figure of a /m n versus m n , we simply drop this logarithmic term which itself vanishes at 
the chiral limit. We also note that for the Wilson fermion action the term oc 1/m 2 may also 
exist, arising from explicit breaking of chiral symmetry, and also from quenching effects [9]. 
We do not see a 1/m 2 effect, as our simulation is perhaps well away from m 2 = and such a 
term is already damped into noise for the range of our simulation. Hence we do not include 
this term into our fit. In order to detect these two additional terms a simulation is needed 
close to the chiral limit with much higher statistics. 

In Fig. 3 we present a /m n in the chiral limit as a function of the lattice spacing, together 
with continuum extrapolations. Their numerical values are tabulated in Table. IV where 
values for (Sqr) are also listed for completeness. This figure demonstrates a sizable scaling 
violation, but exhibits a very clean linear dependence as a function of a. It is interesting 
to observe that the difference between (Exp) and (Lin), which are quite sizable on finite 
lattices, vanishes approaching the continuum limit. This shows that the second order term 
0(t 2 ) included in (3) becomes irrelevant as AE ■ r becomes sufficiently small; one may 
use any formula correct to the first order in r to extract the AE. On the other hand, 
the extrapolation with (Old) gives a value somewhat different from the other two in the 
continuum limit, indicating that the departure of Z from unity could be non-negligible 
(although at 1.2 — 1.5a). 

As our final value for the scattering length in the continuum limit at physical pion mass 
we take the result from (Exp), which agrees with that from (Lin) but has a larger statistical 
error: 



where a rather large error arises from the continuum extrapolation. This result is compared 
with the CHPT prediction : 




2.09(35) 1/GeV 2 , 



(8) 




2.265(51) 1/GeV 2 , 



(9) 
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The scattering length we derived at the continuum limit agrees well with the prediction 
of CHPT. The difference seen in the fitting formula of (Old) and (Lin) accounts for the 
1.5(7 difference of the lattice result from the CHPT prediction mentioned in our preliminary 
report, which is based on an incorrect extrapolation formula (Old). 
We remark that our results also agree with those of Liu et. al. [6] 



where two values (Scheme I and II) refer to their two different treatments for the finite 
volume corrections. 

In this article we have reported a calculation of the scattering length for the I — 2 S- 
wave two pion system. We have shown that the result in the continuum limit is virtually 
independent of the choice of fitting functions used to extract AE from the ratio R(t), and 
that it is consistent with the prediction of CHPT within our 15% statistical error. 



a /m. 
a /m. 



1.75(38) 1/GeV 2 for Scheme I , 
2.34(46) 1/GeV 2 for Scheme II , 



(10) 
(11) 
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FIG. 2. The mass dependence of ao/m n (l/GeV 2 ) at each lattice spacing. 
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FIG. 3. ao/m 7r (l/GeV 2 ) at the chiral limit at each lattice spacing. The CHPT prediction is 
also plotted. 
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TABLES 



(3 = 5.9 










K 


lit 




hj 


/ 

ao/m^ 


2 f/~* \t2\ 

(Lxev ) 




(xlO ) 


/ - , i n — 5 \ 

(xlO ) 


(1/GeV ) 


0.1585 


Old 


12.4(21) 




-0.84(12) 


0.2529(56) 


Exp 


20.9(40) 




— 1.29(20) 




Lin 


14.5(19) 




—0.96(11) 




hqr 


23.1(74) 


on /oi \ 

29(21) 


i /i n / o cr \ 

— 1.40(35) 


0.1580 


Old 


12.5(15) 




-0.822(84) 


0.3468(49) 


Exp 


inn 

19.9(27) 




—1.20(13) 




Lin 


14.0(13) 




—0.905(72) 




hqr 


in n / c t\ 

19.0(57) 


14(15) 


1 1 C* { T7\ 

— 1.16(27) 


0.1575 


Old 


12.1(12) 




—0.786(65) 


0.4396(48) 


Exp 


i o r/oi \ 

18. 5(21) 




i ino/no\ 

—1.108(98) 




Lin 


1 o O / 1 1 \ 

13.3(11) 




—0.849(56) 




hqr 


16.3(50) 


8(13) 


i nn / o /i \ 

— 1.00(24) 


0.1570 


Old 


11.5(10) 




—0.743(55) 


0.5337(49) 


Exp 


i T n / 1 't\ 

17.0(17) 




i n I T / ^rn \ 

—1.017(79) 




Lin 


i o ^ o ( no\ 

12.48(92) 




—0.794(47) 




o , 

bqr 


1 ,1 ^ { A C \ 

14.4(45) 


5(12) 


n on /oo\ 

— O.oy{22) 


0.1565 


Old 


10.86(91) 




-0.698(48) 


0.6297(50) 


Exp 


15.6(15) 




-0.931(67) 




Lin 


11.69(82) 




-0.741(42) 




Sqr 


13.0(41) 


3(10) 


-0.81(20) 


0.1560 


Old 


10.19(82) 




-0.654(43) 


0.7279(51) 


Exp 


14.2(13) 




-0.855(59) 




Lin 


10.92(75) 




-0.692(38) 




Sqr 


11.9(37) 


2.6(95) 


-0.74(19) 



TABLE I. The results at f3 = 5.9. Four lines for each m n are results with the fitting functions 
Old, Exp, Lin, and Sqr, which are defined in (4)-(6). 
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Old 
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JliXp 


11.0(11) 




1 /i / 1 1 \ 
— 1.42(11) 




T in 

Liin 
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1 1 E /I /TO\ 

— 1.154( 1 2) 




oqr 


Iz.z(ol) ) 


1 .Z[oo ) 
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0.153 1 


Old 


/ .38(02 J 




— 0.yb0(b8) 
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Lin 
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oqr 
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0.0(00 ) 


1 oo/oo\ 

— L.ZZ{Z,Z ) 


0.15300 


/■All 

Did 


b.48(4y ) 




—0.831(53) 


0.54b5(4y) 


JliXp 


8.00(08 ) 




— l.OoO(bo) 




Lin 


I .Llyil ) 




— u.oyo(4o ) 




Sqr 


y.o(zi j 


4. < (4y j 


1 1 1 f 
— 1.11(ZU ) 


0.15250 


/All 

(Jlu 


5.yb(43 J 




— 0. r 00(45 ) 


0.b8y< (52) 


JliXp 


1 . 1 0(08 ) 




— 0.y38(5b ) 




Lin 


6.48(42) 




-0.814(42) 




Sqr 


8.2(19) 


3.9(43) 


-0.99(18) 


0.15200 


Old 


5.48(40) 




-0.697(41) 


0.8385(55) 


Exp 


6.95(52) 




-0.845(50) 




Lin 


5.92(38) 




-0.743(39) 




Sqr 


7.4(17) 


3.1(39) 


-0.89(16) 



TABLE II. The results at f3 = 6.1. Four lines for each m n are results with the fitting functions 
Old, Exp, Lin, and Sqr, which are defined in (4)-(6). 
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Sqr 
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0.8(19) 
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Old 


4 36(30") 




-0 842(48") 


0.6634(45) 


Exp 


5.37(39) 




-0.996(58) 




Lin 


4.70(30) 




-0.894(46) 




Sqr 


4.72(89) 


0.0(17) 


-0.90(14) 


TABLE III. 


The results at f3 


= 6.3. Four lines for each 


are results with the fitt 


ing functions 


Old, Exp, Lin, and Sqr, which 


are defined in (4)-(6). 
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(3 


a (1/GeV) 


Old 


Exp 


Lin 


Sqr 


5.9 


0.493(7) 


-0.96(10) 


-1.51(16) 


-1.093(90) 


-1.58(36) 


6.1 


0.378(6) 


-1.185(59) 


-1.653(80) 


-1.335(55) 


-1.78(22) 


6.3 


0.302(5) 


-1.335(76) 


-1.745(99) 


-1.466(74) 


-1.77(21) 




a -► 


-1.92(25) 


-2.09(35) 


-2.07(24) 


-2.04(78) 



TABLE IV. The values of ao/m^l/GeV 2 ) in the chiral limit for the each fitting function of 
R(t) at each (3 and those in the continuum limit obtained by liner extrapolation in the lattice 
spacing. The fitting functions of R(t) are defined in (4)-(6). 
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